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Dissipative Collapse of a Spherical Cluster of Gas Clouds 
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Abstract. With a view to understand the galaxy /star formation scenario, we investigate the dissipative collapse 
of a spherical cluster of gas clouds with an isotropic velocity distribution. The time scale for collapse to one tenth 
radius is studied as a function of the collision time in the system. The scalar virial equation is used to investigate 
the evolution of the size of the cluster. This is supplemented with an evolution equation for the random kinetic 
energy. The above system is numerically solved and the results analyzed. For large values of the collision time we 
find that the time scale for collapse is proportional to the collision time as expected. However for large values of 
the dissipation, i.e. for small collision times, the collapse time shows a nonlinear dependence on the collision time. 
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1. Introduction 

The formation of stars and galaxies normally proceeds 
from the collapse of gas clouds. Hence it is imperative that 
one investigates the evolution of the size and eccentricity 
of spherical and spheroidal clusters of gas clouds. Such an 
investigation for point masses with an isotropic distribu- 
tionof velocities was made by Chandrasekhar and Elbert 
(1972) using the scalar and tensor virial equations. This 
was followed by Som Sunder and Kochhar (1985), who fur- 
ther examined (1986) the evolution for anisotropic velocity 
distributions also. Here we examine the dissipative evolu- 
tion of a spherical cluster of gas clouds with an isotropic 
velocity distribution, using the scalar virial equation, in 
an attempt to obtain the time scale, for collapse impeded 
by virialized mass motions. The gravitational binding en- 
ergy released during collapse feeds the random motions 
of the clouds. This kinetic energy is dissipated in cloud- 
cloud collisions. The evolution of the cluster size is fol- 
lowed using the scalar virial equation. It is supplemented 
with an evolution equation for the kinetic energy of ran- 
dom motions. From this investigation we get quantitative 
relations between the collapse time and the collision time 
under various conditions. 



formly. The clouds have a one dimensional r.m.s. velocity 
Vrms- The gas clouds are in pressure balance with an in- 
tercloud medium. The mass in the intercloud medium is 
taken to be too small to affect the dynamics and the con- 
tribution from it to the virial as well as the dissipation 
will be ignored throughout (Mestel and Paris 19841 ). The 
scalar virial equation for the system is (Chandrasekhar 
and Elbert |1972| ) 



2T + n 



(1) 



2 dt^ 

where I — Att J p{r)r^dr is like the moment of inertia 
of the system about the centre and p{r) is the density at 
a distance r from the origin. T is the kinetic energy asso- 
ciated with mass motions, and is the self gravitational 
potential energy of the cluster. The evolution of the sys- 
tem must be consistent with equation 1 and we may use it 
to investigate the evolution of the gross size of the cluster. 
For the cluster we get 

/=^Mi?^(l + (f )^) (2) 
and 



2. The model 

We consider a spherically symmetric cluster of A'' equal 
mass clouds of radius Rr and mass Mr distributed uni- 
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(3) 



Thus for Rc much smaller than R using the expressions 
for I and fl corresponding to a homogeneous distribution 
of matter would not be too erroneous. 
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The velocity of the clouds may be separated into two 
parts, a random part and a mean motion. The kinetic 
energy for the random motions say Tr is |Mw^„j^. If we 
consider only homologous changes in the size of the system 
with the mean velocity at a point being proportional to the 
distance of the point from the centre, we get the kinetic 
energy for the mean motion as Tm = Jo^^i^)"^ ■ The 
cluster is assumed to be uniform at all times. Introducing 



and 



H f) "rms 



the virial equation becomes 
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The evolution equation for the random kinetic energy 
will have a growth term due to the feeding of gravitational 
potential energy into the random motions and a dissipa- 
tive term due to cloud collisions. In the case of nonuniform 
collapse the gravitational energy released during contrac- 
tion goes initially into radial flow, but eventually some of it 
is converted into the energy of random motions. The possi- 
bility that all the released binding energy is converted into 
random kinetic energy has also been advocated (Odgers 
and Stewart 1958 ). We consider that all the released po- 
tential energy is converted to random kinetic energy. This 
gives 



dTR _ dn 

\growth~ 



dt 



dt 



(7) 



For dissipation by collisions the dissipation rate will 
be proportional to (the number of collisions suffered by 
a cloud in unit time) x (the energy loss per collision) x 
(the total number of clouds). The number of collisions per 
unit time will be proportional to the (geometrical cross 
section of the cloud) x (the velocity of the cloud) x (the 
number density of clouds). The dissipation rate is thus 
proportional to Mv^^^NRl/ R^ =vMvl^^/ R. Here we 
have expressed Rc in terms of the cloud filling factor f 
as Rc = {jf)3R. Putting all the numerical factors into a 
parameter a we write i/ = aN^ fs . Here v is proportional 
to the number of collisions in a free fall time for virial 
theorem random motions in the system. As an example. 



Elmegreen(1985) when considering that the initial speed 
of a cloud is halved after interacting with a mass equal to 
itself via collisions, gets a = We thus get 



Putting all the terms together we get the evolution 
equation for the random kinetic energy Tr in the form 



dvt: 



dt 



-2GM dR 
5i?2 'dt 



2v^ 
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(9) 



As the cluster contracts it is the contraction produced 
by the kinematics of the collapse itself that mainly af- 
fects the radius of the clumps rather than any self grav- 
itational effects Lacey(1984). For clouds which are them- 
selves Jeans' unstable the collapse time for the clouds will 
be smaller than the free fall collapse time for the cluster 
only by a factor f. The number of the clouds can vary as 
there can be cloud disruptions or coalescences also during 
collisions. We assume a constant v throughout the evolu- 
tion. 

For corresponding to virial equilibrium the 

collision time is ~l/i/ in units of the free fall time 

^ff^VioM)^ where i?o is the initial radius. Thus since 
dissipation is due to collisions the dissipation time is also 
^l/v. For the system to be coUisional the mean free path 
of the clouds should be less than twice the radius of the 
cloud. This yields the condition Ns fa > 1/6. This condi- 
tion is independent of the size of the cluster since we have 
assumed a constant filling factor throughout the collapse. 



3. Results 

We have considered evolution of the cluster under two 
different conditions: (a) the system is initially in virial 
equilibrium with Tji = fl/2 (b) the system has zero ini- 
tial kinetic energy. As case (c) we have considered the 
evolution of the cluster starting from virial equilibrium 
and with no gravitational feeding . Equations (6) and (9) 
were solved for various values of the only free parameter 
v^l/tc where tc is the collision time for virial theorem 
random motions in units of tff. The equations were nor- 
malized and integrated with a constant time step. Desired 
accuracy was achieved by using a time step which was 
appropriately small. It was checked that with zero dissi- 
pation the change in the total energy was less than one 
percent over 10 units of normalized time. 

In Fig. 1 we give R as a function of time for the cases 
(a) and (b) for two different values of i^, .2 and 0.5. The 
two curves with faint oscillations are for case (a). The 
curves with strongly marked oscillations are for case (b) 
and the curves for no gravitational feeding, case (c) show 
monotonic collapse. For larger values of i' the oscillations 
start at smaller radii. 

In Fig. 2 we show tio the time for collapse to one-tenth 
size in units of tff as a function of tc{— {aN^^^ 
The uppermost curve is for case (a), the lowermost for 
case (b), and the middle curve is for case (c). For large 
values of ic, ^lo depends linearly on tc and a fit to the 
curves in this region gives, for case (a) 



dTR 
dt 



\loss ~ 



R 



(8) tio^0.95 + 2.0tc 



(10) 
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Fig. 1. R as a function of t for tc=5 and 2 for cases(a), 
(b) and (c) with t in units of tff. 
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Fig. 2. as a function of tc for cases (a),(b) and (c) 
Both times are in units of tff. 

and for case (b) we get 



0.74 + 0.6tc 



(11) 



In the nonlinear region the curves can be represented 
with very good accuracy by a quadratic in tc- For case (a) 
we get 



tio = 1.25 + 3.0tc- lAtl 
and for case (b) 



1.16 + 0Mtc-0A6K 



(12) 



(13) 



Fig. 2 also shows the case of no gravitational feeding 
(case (c)), for comparison with earlier work where the dis- 
sipation of virial theorem random motions was investi- 
gated by various authors. The cluster is initially in virial 
equilibrium. For small dissipation i.e. for large collision 



times tio is not very sensitive to tc and increases slowly 
with tc in a linear fashion as tio — 2.2 -I- 0.2ic- However 
for large dissipation, there is a fast increase of the collapse 
time as the collision time increases. This nonlinear portion 
is fitted by the curve t^ = 1.2 + 2.1tc - 1.38*^. 

4. Discussion 

For a nondissipative system starting with zero kinetic en- 
ergy Chandrasekhar and Elbert ( 1972| ) have shown that 
the system will execute oscillations. If we put zero dissi- 
pation in our system we get the same oscillations. For the 
finite dissipation cases which we have considered we find 
these oscillations; however they are of a damped character 
due to the finite value of the dissipation. 

Examining i? as a function of t, in all the cases with 
gravitational feeding we find that collapse takes place in an 
oscillatory fashion. The oscillations are well marked in case 
(b) . Here there is an initial rapid collapse when Tji rapidly 
builds up and finally reaches a value high enough for the 
virial equilibrium condition to hold. This takes place when 
the value of ^ is quite large and the resulting bounce of 
the system is quite hard. This causes the oscillations to be 
strongly marked. As v increases the bounce is delayed due 
to the the slower increase of Tr. In case (a) the bounce 
can take place quite early on and is quite mild. For large 
dissipation the bounce gets delayed so much that we do 
not get it within collapse to one tenth size and the collapse 
looks apparently monotonic. 

It is known that a self gravitating cloud collapsing 
without dissipation will satisfy the conditions for virial 
equil ibrium when it has collapsed to half radius (Hoyle 
1949 ). The idea of virialization of mass motions in gravi- 
tational collapse and their subsequent dissipation has been 
conside red b y many people. (Odgers and Stewart 
Larson 1977, M estel 1977 , Scalo and Pumphrey 
Elmegreen 1985.) These works consider virial theorem 
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1982 



random motions and the time scale for their dissipation. 
Simulations of the coUisional dissipation of virial theorem 
random motions by Scalo and Pumphrey give a dissipa- 
tion time ~20 f ree fa ll times. (See also Chieze and Lazareff 
|1980| , Hausman [l98l| .) Elme green considered a magnetized 
cloud of clumps and obtained dissipation times of the or- 
der 3 — 5 cloud crossing times. 

We have analyzed collapse with gravitational feed- 
ing, of clusters starting from virial equilibrium and also 
starting with zero kinetic energy (like for example from 
turn around) . For small values of the dissipation like say 
tc — 5tff we get tio about 11 for case (a) and about 3.8 
for case (b). In the model of Scalo and Pumphrey (1982) 
the number of clouds get smaller with time due to cloud 
coalescence. The resultant slower dissipation rate explains 
the comparatively longer dissipation times they get. As a 
general statement we may say that for small dissipation, 
clusters starting from virial equilibrium can take more 
than twice the dissipation time to collapse, when there 
is gravitational feeding. With no gravitational feeding we 
get > tff > tc- The increase in the collapse time as 
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the collision time increases is less pronounced in this case Som Sunder G., Kochhar R.K. 1985, MNRAS, 213, 381 
and for large collision times lies in the range 2-4 free Som Sunder G., Kochhar R.K. 1986, MNRAS, 221, 553 
fall times. This is consistent with the results obtained by 
other authors for the dissipation time for virial theorem 
random motions. As tc goes to zero tio rapidly approaches 
the value 1.1 which is the time for free fall collapse to one 
tenth size in the units we have chosen. Our results in- 
dicate that for large collision times gravitational feeding 
can significantly increase the time scales for the dissipative 
collapse of clusters of gas clouds. 

For large values of tc i.e. small values of the dissipation 
we sec from fig. 1 that the collapse time depends linearly 
on tc- For case (a) this may be understood as follows. 
For small dissipation the cluster is collapsing very slowly. 
Thus growth as well as loss of Tr is small and the clus- 
ter contraction is governed mainly by the decrease of Tr 
in the virial equation. This decrease is in proportion to 
tc- For case (b) also after an initial rapid collapse, when 
Tr rapidly builds up, the size of the cluster is governed 
similarly by the decrease of Tr in equation (6). For large 
dissipation however the collapse being rapid, changes in 
^ and R are large causing large changes in Tr. There is 
thus a strong and nonlinear dependence of the evolution 
of Tr on tc- The collapse time also reflects this strong 
nonlinear dependence. 

For the system under consideration the free fall time 
provides the appropriate time scale for the system to ad- 
just to changes in the virial theorem (Brosche 1970). While 
noticing that a quasi virial evolution does not hold for dis- 
sipation times shorter than the dynamical time, we find 
distinctive difference in the dcpcndance of tio on tc in the 
two regions tc > tfj and tc < tfj. 

The consequences of this investigation for the forma- 
tion of stars and galaxies are being examined at present 
in a quantitative fashion and will be reported separately. 
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